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ABSTRACT 

We analyse a sample of 507 evolved (OH/IR) stars in the region (10° > £ > 
— 45°) , < 3°). We derive average ages for subsets of this sample and use those 
sets as beacons for the evolution of the Galaxy. In the Bulge the oldest OH/IR stars 
in the plane are 7.5 Gyr (1.3 M Q ), in the Disk 2.7 Gyr (2.3 M ). The vertical dis- 
tribution of almost all AGB stars in the Disk is found to be nearly exponential, with 
scaleheight increasing from 100 pc for stars of ^1 Gyr to 500 pc for stars of ^5 Gyr. 
There may be a small, disjunct population of OH/IR stars. The radial distribution of 
AGB stars is dictated by the metallicity gradient. Unequivocal morphological evidence 
is presented for the existence of a central Bar, but parameters can be constrained only 
for a given spatial-density model. Using a variety of indicators, we identify the radii of 
the inner ultra-harmonic (2.5 kpc) and corotation resonance (3.5 kpc). We show that 
the 3-kpc arm is likely to be an inner ring, as observed in other barred galaxies, by 
identifying a group of evolved stars that is connected to the 3-kpc HI filament. Also, 
using several observed features, we argue that an inner-Lindblad resonance exists, at 
~1-1.5 kpc. The compositions of OH/IR populations within 1 kpc from the galactic 
Centre give insight into the bar-driven evolution of the inner regions. We suggest 
that the Bar is ~8 Gyr old, relatively weak (SAB) and may be in a final stage of its 
existence. 
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1 INTRODUCTION 

It has become widely accepted that our Galaxy is barred, 
as evidence accumulated over the last five years from star 
counts, gas-dynamical studies and stellar three-dimensional 
kinematics and especially from the analysis of the COBE- 
DIRBE integrated-light data (Dwek et al. 1995; for a re- 
view on the galactic Bar see Gerhard 1996). The important 
parameters of a barred potential are the semi-major-axis 
length, a, the in-plane axis ratio, q, the pattern speed, fi p , 
and the strength relative to the axisymmetric part of the po- 
tential, A. From the observer's point of view, another quan- 
tity is the major axis' orientation with respect to the line 
of sight, the viewing angle <j>- Due to the awkward view we 
have of the Galaxy, the values of those parameters are even 
harder to determine than in external galaxies. They are of- 
ten inferred from the influence the barred potential has on 
the other parts of the galaxy. Helpful, albeit crude and not 
completely understood, diagnostics are the "resonant rings" 
and spiral features that arise at radii where orbits are in res- 
onance with the frequency of rotation of the bar. The use of 
these structures is especially difficult in our Galaxy, where 
we see only tangent points to rings and spiral arms, but we 



will attempt to do so. 

The aim of this article is to describe the overall form of 
the stellar distribution in, and the evolution of various galac- 
tic components (disk, bulge, spiral arms etc), based mainly 
on its content of evolved stars and focussing mainly on the 
inner Galaxy. To a lesser extent, we attempt to constrain 
the free parameters of these stellar distributions. The goal 
is a schematic, rather than comprehensive, picture of the 
Galaxy. 

We restrict ourselves largely to the stellar distribution 
and discuss the distribution of the gas only superficially. As 
stars are probably the source and the cleanest tracers of the 
barred potential, such a study of only the stellar compo- 
nent of the Galaxy is necessary. Obviously, the gas- and the 
stellar distributions should ultimately be explained simulta- 
neously in one coherent picture. In fact, the first evidence 
for non-axisymmetry of the inner Galaxy came from the 
neutral-gas kinematics (see review by Oort 1977). It was 
exactly the drive for coherence that has shifted the atten- 
tion to the stars, because their distribution showed no clear 
deviation from axisymmetry. Neither of the prevailing ex- 
planations for the origin of the observed radial gas motions 
- central expansion or elongation of the potential - were sup- 
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ported by stellar observations. Once the Bar had been found 
in the stellar surface density (Blitz & Spergel 1991; Dwek 
et al. 1995), questions remained. For instance, why are the 
stellar kinematics explained so well by axisymmetric models 
(Kent 1992; Ibata & Gilmore 1995) and why is the micro- 
lensing optical depth toward Baade's window incompatible 
with density models derived from surface-density maps (eg. 
Nikolaev & Weinberg 1997) ? 

The stellar data used in this paper (Sevenster et 
al. 1997a,b) were collected specifically to form the optimal 
sample to address such issues and to complement existing 
data. The sample consists of OH/IR stars: intermediate- 
mass, oxygen-rich, far-evolved asymptotic-giant-branch 
(AGB) stars. These are excellent tracers of the general stel- 
lar population, as stars with initial masses between roughly 
1 M and 6 Mq go through this phase. They are also tracers 
of the galactic potential as they form a fairly relaxed popula- 
tion, with typical ages of several gigayears. Very strong and 
characteristic maser emission at 1612 MHz from the ground- 
state OH molecule allows for radio-interferometric observa- 
tions, combining extinction-free coverage of the plane and 
fast sampling out to large distances with positions and veloc- 
ities with negligible errors (compared to modelling errors). 
The employability of this sample in galactic-structure stud- 
ies is clearly demonstrated in Sevenster et al. 1999. 

The structure of this paper is as follows. In the first part, 
we discuss the large-scale spatial distribution of the OH/IR- 
star sample. We describe the sample of OH/IR stars in more 
detail in Sect. 2. In Sect. 3 we analyse the structure of this 
sample and estimate the scaleheights and -lengths for the 
galactic Bulge and galactic Disk. Variations of those expo- 
nential scales with position are combined with scales derived 
from similar samples and interpreted as age dependencies. 
In Sect. 4, we give morphological evidence for a galactic Bar 
that cannot be explained by a physical lopsided distribution 
(see Blitz & Spergel 1991; Sevenster 1996). By comparing 
some parametrized models to the observations, a flavour of 
the values of the relevant parameters is obtained. Most im- 
portantly, we derive an approximate distance limit for the 
sample. The kinematic type of the galactic Bulge in terms 
of its anisotropy parameter and ellipticity is found. 

In the second part of the paper we discuss the influ- 
ence of the bar-shaped potential on the inner Galaxy. In 
Sect. 5 we interpret various patterns found in the sample as 
resonant patterns and identify resonant radii. New insight 
into the nature of the so-called 3-kpc arm is provided by a 
small group of OH/IR stars that follows the 3-kpc arm's 
longitude-latitude-velocity structure. In Sect. 6 we give a 
description of the inner Galaxy based on the presented evi- 
dence and we speculate on its evolution in Sect. 7. We con- 
clude in Sect. 8. 

Throughout this article, we will use the term "Bulge" 
to denote the galactic component we see in the general di- 
rection of \£\ £,10°, without being interested in its actual 
form. If we use the term "Bar" we denote specifically the 
prolate or triaxial component of the Bulge. Unless explicitly 
stated otherwise, stellar velocities are given with respect to 
the local standard of rest (see Sevenster et al. 1997a,b for 
the used Doppler corrections) and Rq =8 kpc. 



2 PROPERTIES OF THE DATA SET 

The stellar density and the gravitational potential of the 
Galaxy are best traced by intermediate-mass, evolved stars 
as they constitute the largest fraction of the total stellar 
mass and are dynamically relaxed (Frogel 1988). Good can- 
didates for this are the so-called OH/IR stars; oxygen-rich 
objects on the asymptotic giant branch (AGB; see for in- 
stance Habing 1993; Sevenster, Dejonghe & Habing 1995). 
We use an unbiased, homogeneous sample of 507 OH/IR 
stars in the galactic Plane, the AOSP sample (Australia tele- 
scope Ohir Survey of the Plane), acquired in a systematic 
survey of the region between longitudes 10° > I > —45° 
and latitudes |6| < 3° in the 1612 MHz (18 cm) OH-maser 
line (Sevenster et al. 1997a,b). The positional accuracy is 
0'.'5, the line-of-sight velocities (with respect to the local 
standard of rest) are accurate to 1 kms -1 . 

In Fig. 1(a) we show the surface density of the AOSP 
sample, smoothed with an adaptive-kernel algorithm (Mer- 
ritt & Tremblay 1994). With the initial kernel size that re- 
tains best the steepness of the central density profile with- 
out showing individual stars (1°) many local maxima are 
revealed in the distribution (R-i, R-2 and R_3). The max- 
imum at t = +9° is likely to be a spurious edge effect. We 
will not discuss any asymmetries in the vertical direction, 
such as a tilt, because the data quality is slightly latitude- 
dependent (Sevenster et al. 1997a, b). First, in this section, 
we will assess some astrophysical properties of the sample. 

The OH emission comes from an optically-thick, ex- 
panding circum-stellar envelope. Because the expansion is 
radiation-pressure driven, the unobservable intrinsic stel- 
lar luminosity L* is related to the outflow velocity of the 
circum-stellar envelope, Kx P , and its gas-to-dust ratio, 
(related to metal abundance Z\ \i oc Z -1 for oxygen-rich 
stars, Habing, Tignon & Tielens 1994). According to van 
der Veen (1989): 

L, oc ^i 2 Kt p . (1) 

This equation is derived and discussed in more detail in Ap- 
pendix A. Separating the objects according to outflow veloc- 
ity (Fig.l(b,c)) hence results in a separation roughly accord- 
ing to stellar luminosity, or age, without knowledge of the 
distances to the objects. This results, as expected, in very 
different apparent scaleheights of the subsamples (Fig.2, see 
discussion in Sect. 3). 

OH/IR stars can span a wide range of ages, of 0.1 Gyr 
to <;10 Gyr. From the IRAS two-colour diagram, one can 
determine a "turn-over" [25]-[60] colour R ! 32 , of sources 
leaving the evolutionary track (van der Veen & Habing 
1988). This R 32 is related to the initial mass (Garcia Lario 
1991); this relation is discussed in Appendix B. We can thus 
find the initial mass of a star that has reached the end of 
the OH/IR-star phase (the thermally-pulsing AGB phase). 
Since this phase itself is short (~ 10 5 yr, Tanabe et al. 1997) 
compared to the ages of the AGB stars, we can use the 
value for the AGB-tip ages by Bertelli et al. (1994) to ob- 
tain, for an assumed abundance, the age of the star from 
its initial mass. The minimum R 32 for a sample hence gives 
the maximum age of the stars in that sample. We deter- 
mined this for the Bulge region {R 32 > —0.4, Sevenster et 
al. 1997a,b) and the Disk region (R 32 > +0.2, Sevenster et 
al. 1997a,b) separately. We find minimum initial masses and 
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Longitude 

Figure 1 The surface density of the AOSP sample (Sect. 2), in longitude versus latitude in degrees. The observations were smoothed with 
an adaptive-kernel algorithm (initial kernel of 1°). The upper panel is for all OH/IR stars (507), the middle for outflow velocities higher 
than 14 kms" 1 and the lowest panel for outflow velocities lower than 14 kms -1 (and higher than 1 kms~ 1 ). Panels b, c represent 
younger and older stars, respectively, with an abundance effect (Sect. 2). Both subset-plots are based on ~210 stars. Note the offset 
toward negative longitudes in for the young sample (see Sect. 4.1). The local maxima at ~ —10°, —20°, —30° - R_ i, R_2 and R_ 3, 
respectively — will be discussed in Sect. 5. The unmarked vertical dashed lines indicate maxima in the observed 2.4 GHz distribution (see 
Sect. 5, Fig.ll). The maximum at £ = +10° is not reliable. Contours are spaced at twenty even intervals between zero and the maximum 
of each particular plot. 
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Latitude 

Figure 2 The cumulative-number distribution over latitude for 
the two subsamples of the AOSP sample (Fig.l(b,c)). The (ap- 
parent) scaleheights differ by a factor of 2, which is indicative of 
the age difference between the samples. 



maximum ages of 1.3 M©& 7.5 Gyr and 2.3 M©& 2.7 Gyr, 
respectively, both for solar abundance. The ratio of the mean 
initial masses (3.3/2.3 with an IMF oc M~ 2,5 , almost inde- 
pendently of Z used to derive the ages) and of mean outflow 
velocities (15/14 kms -1 ) are compatible with a similar av- 
erage abundance for Disk and Bulge stars, Zb = 0.9-0.95Zd, 
according to equation (1). The number of stars in the AOSP 
sample with ages SS0.5 Gyr is negligible. 

These values are meant to give an indication only, with 
a likely uncertainty of about 15% in the derived ages and 
about 10% in the initial masses. The R^ 2 — Mi relation is not 
yet fully established and should be used only for ensembles, 
not for individual objects. Besides, not all AOSP stars have 
a (reliable) IRAS identification (Sevenster et al. 1997a,b), 
although those that do must be representative of those that 
don't, as the identifications are hampered mainly by the 
confusion-limited spatial resolution of IRAS. 

The ages are influenced by the fact that the AOSP sam- 
ple covers only low latitudes. The fraction of (sub-) solar- 
mass OH/IR stars is apparently smaller than 1% at low 
latitudes. Out of the plane, van der Veen & Habing (1990) 
find the masses of oxygen-rich AGB stars to range to < 1.0 
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Mq, and the corresponding ages to well over 10 Gyr. How- 
ever, with their mass-loss parameter =1.5 instead of 
=1.0 (both values are equally likely to be right), the mass 
range would be 1.2-2.2 Mq, which fits in very well with our 
derivations. Correspondingly, the ages of their stars would 
be lower, especially when using Z = 0.02 instead of Z = 0.04 
as they advocate. We conclude that our results are compat- 
ible with theirs. 



3 EXPONENTIAL SCALES OF THE DISK 
AND THE BULGE 

The latitude distribution of the AOSP sample is steeper 
than a (projected) exponential (exp(— z)) at all longitudes 
and irreconcilable with flatter functional forms (sech 2 (z) 
or exp(-z 2 )). The same was found recently by de Grijs & 
Peletier (1997) for a large number of spiral galaxies; Kent, 
Dame & Fazio (1991) found from NIR observations of the 
Galaxy that the profile is closer to exponential than to 
sech 2 (z). To obtain estimates of the values of the scalelength 
and scaleheight of the spatial density, we hence used double 
exponentials (p oc exp( — R/h-R.) exp(— z/h z )). We estimated 
those values in various regions of the sky, to assess variations 
of the scales with radius and height (age). Note, however, 
that we are not trying to model the Galaxy as a set of double 
exponentials. 

We fitted cumulative-number densities with single ex- 
ponentials, in latitude (in R tan6); Rq = 8 kpc) and lon- 
gitude (in R sin £) separately. This yielded apparent scales 
hzb and h^t of the local surface-density distribution at vari- 
ous longitudes and latitudes (Fig. 3). Due to the non-trivial 
angles between lines of sight, the scales are not invariant 
for projection, even though the Galaxy is seen edge-on. 
This means we have to deproject the apparent scales to find 
the true scales h z and hn- By projecting analytic double- 
exponential distributions with a range of hn and h z , we ob- 
tained a range of apparent scales at the same (£, b) as the 
data measurements. We could thus retrieve the intrinsic (h z , 
/ir) that would yield a measured pair (h zb , h-g,t). Projection 
effects cause the relation between the true and the apparent 
scales to depend upon longitude and latitude, and also upon 
each other, thus jeopardizing a unique deprojection. 

However, for all our measurements, except those at C 
(Fig. 3), this technique gave unique results (Table 1, columns 
3-6), within the context of using double exponentials. In Ta- 
ble 1 the deprojected scales are also given for the high- and 
low-outflow sources separately, determined over the whole 
latitude range of the survey. For the galactic disk (|^|^15°), 
most references (see Sackett 1997) give hn =2.5-4.5 kpc 
and h z =250-400 pc ("thin") and h z =750-1500 pc ("thick" 
disk). 

3.1 Ring ? 

The very large apparent scalelength at I ~ —18° (C; Fig. 3) 
indicates a flat distribution in longitude. With this value of 
the measurement, the deprojection is not unique. Neverthe- 
less, our tests showed that projection effects alone, although 
being largest at these longitudes, can not explain the ex- 
treme apparent scalelength and we conclude it is intrinsic. 
The flat distribution is formed mainly by the high-outflow 
(and single-peaked) sources, ie. the younger stars (Table 1, 
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Figure 3 The local apparent exponential scaleheights h zb and 
h-ni of the AOSP sample, from fits to the cumulative distribu- 
tions in latitude and longitude separately. The symbols are for 
< |&| < 1 (closed) and 1 < |6| < 3 (open), respectively. The 
bars along the abscissa indicate the longitude bins (A,B,C,D) 
used for the determination of the fits; the data points are at the 
middle values of these bins. All longitude bins contain about 100 
stars over the full latitude range (\b\ < 3°); the formal errors on 
the measurements arc ~15%. The two circles are the scalclcngths 
fitted to low-outflow stars (cf. Fig. 1(c)) only. For A the apparent 
scaleheight is very small due to the presence of the galactic-centre 
population (see Sect. 3. 2). Very large apparent scalclcngths, such 
as in longitude bin C, are equivalent to an infinite scalelength or 
a (locally) flat distribution (Sect. 3.1). 

column 7 vs. 9; Fig. 3 (circles)). This region around I = —18° 
will be treated in more detail in Sect. 5, where we find that 
the distribution in that direction is probably dominated by 
a ring structure. 

3.2 Central disk 

A deciparsec-scale, flat, rapidly-rotating flat population of 
OH/IR stars is present in the galactic Centre (Lindqvist, 
Habing & Winnberg 1992; Sevenster et al. 1995). Even 
though the AOSP sample contains only 19 out of the 134 
known OH/IR stars in this region (Lindqvist et al. 1992; 52 
more were discovered by Sjouwerman et al. 1998a), we can 
still resolve its very low scaleheight with respect to other re- 
gions (Table 1, column 5). As was known before, the central 
disk is seen primarily in the high-outflow sources (Table 1, 
column 10). 

However, the scalelength of 250 pc (Table 1, column 3) 
is unlikely to represent this small central disk. Indeed, if we 
apply our method to the full sample of Lindqvist et al. (1992) 
we find h z = 20 pc and /ir = 40 pc, the same values as they 
find themselves (Table 1, columns 3,5 between brackets). 



3.3 Outer Galaxy 

Hardly any OH/IR stars are known outside the solar cir- 
cle. Carbon-rich AGB stars, on the other hand, are hardly 
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Table 1. Deprojccted exponential scales for the squares and triangles in Fig. 3 (columns 3—6) and for low- ("o", 
columns 7,8) and high-outflow ("y", columns 9,10) sources (|6| < 3°). The values between brackets are from another 
sample at \£\ = 0° (see text) or conceivably unreliable because of the infinite value of at \£\ = 18°. 
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Figure 4 The cumulative-number distributions of outflow ve- 
locities for double-peaked OH/IR stars from various samples: 1. 
Disk AOSP sample; 2. Arecibo sample (Chengalur et al. 1993); 3. 
Outer Galaxy (Blommaert et al. 1993); B. Bulge AOSP sample; 
GC. galactic— centre sample (Sjouwerman ct al. 1998a). The me- 
dian outflow velocity correlates with metallicity (equation (1)); 
its decrease with larger mean radius (in order : (GC,B) 1,2,3) 
reflects mainly the galactic metallicity gradient. The small dif- 
ference between "1" and "B" is mainly due to a mean-mass dif- 
ference, however. The arrows indicate deviant sub-populations 
of distribution "2" (thick disk; discussed in Scct.3.4) and "GC" 
(star burst; Sect. 7). Except sample "3" that consists of 9 stars, 
all samples contain more than 200 stars. 

seen in the inner Galaxy (see the figures in Blanco 1965). 
Blommaert, van der Veen & Habing (1993) find that the few 
OH/IR stars in the outer Galaxy have initial masses of at 
least 2-3 M , similar to the AOSP disk stars (Sect. 2). The 
low outflow velocities of the sample (the median is 3 km/s 
lower than for the Disk AOSP sample, Fig.4) indicate low 
metallicity (equation (1)), as was observed by Blommaert et 
al. (1993). At intermediate longitudes (£ ~ 50°), the OH/IR 
stars from the Arecibo survey (Chengalur et al. 1993) have 
initial masses and ages similar to the Bulge AOSP sample 
(Sect. 2) and slightly lower median outflow velocity. 

The decrease in outflow velocities with higher longi- 
tude is the result of decreasing metallicity. The galactic-disk 
metallicity gradient (—0.07 dex kpc -1 from oxygen, Smartt 
& Rolleston 1997) corresponds, for constant mass, to an 
outflow-velocity gradient of ~1 kms -1 kpc -1 (for t4 xp ~14 
kms -1 , equation (1)). 

At the same time, at lower metallicity, the limiting 
mass above which stars remain oxygen-rich and below which 
carbon-rich stars form increases. The upper limit to the 



mass for which a star will reach the AGB in turn becomes 
lower with lower abundance (see Bertelli et al. 1994). The 
upper limit to the outflow velocities will decrease even more 
noticeably (equation (1)); this explains the small scalelength 
for the high-outflow sources (Table 1). The abundance may 
decrease rather suddenly around R = 6 kpc, as observed by 
Simpson et al. (1995). They explain this by assuming this is 
the outer edge of the zone-of-influence of the Bar. As we 
will see later (Sect. 6), the Bar's outer-Lindblad resonance 
(OLR, Binncy & Tremaine 1987 (BT) Ch.6) may indeed 
have a radius of 6-7 kpc. 



3.4 The Arecibo sample 

With respect to an otherwise equivalent sample (te Lintel 
et al. 1991), the Arecibo sample has twice as many sources 
at outflow velocities Ss9 kms -1 (Fig.4, see left-most ar- 
row ). This difference is significant to 99.8% (Kolmogorov- 
Smirnov). The apparent scales of the subset with outflow 
velocities < 10 kms~ (Fig. 6) are h z b = 1-1 kpc and 
h Rt — > oo (Fig. 5). Again, the nearly-infinite apparent scale- 
length makes unique deprojection impossible, but we found 
that the apparent scaleheight is most probably close to the 
true scaleheight in this case. These stars might be part of 
the 1-kpc thick disk claimed by other authors (Habing 1988; 
Gilmore & Reid 1983, Ojhaet al. 1996), although Blommaert 
et al. (1993) conclude that true AGB stars do not partake 
in the thick disk and there is no sign of this population in 
the te Lintel (et al. 1991) sample. In Fig. 6, we see that these 
low-outflow stars are not concentrated toward the plane. 
This may be a selection effect as the Arecibo sample was 
acquired in a detection experiment toward IRAS-selected 
point sources, but this is also true for the te Lintel (et 
al. 1991) sample. Moreover, when we select stars from the 
Arecibo sample in the same latitude range as the AOSP 
sample a small excess is still present (Fig. 7). 

So, in the Arecibo sample we see an extra population 
with low outflow velocities and a gap in the latitude distri- 
bution. The latter brings to mind the intriguing "levitation" 
process (Sridhar & Touma 1996a,b), a vertical-to-radial res- 
onance. This would not explain an excess of sources, how- 
ever. If the excess sources have low metallicity, rather than 
only low masses (equation (1)), this population could have 
been accreted anytime during the last 10 Gyr from a dwarf 
such as Fornax ([Fe/H] < —0.7). With the observational 
limit of less than five such more metal-rich dwarfs accreted 
in the last 10 Gyr (Unavane, Wyse & Gilmore 1996) this is 
a possibility. 
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Figure 5 The cumulative-latitude distribution of all Arecibo stars with outflow velocity <10 kms — 1 (solid, cf. Fig. 4, 6). The fit (dashed) 
is a projected double-exponential distribution with apparent scaleheight =1.1 kpc. For comparison, the dot-dashed lines give the 
distribution for h zi) =800 pc and 1 kpc, respectively. The intrinsic scaleheight is probably close to 1.1 kpc as well. 
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Figure 6 The distribution of the Arecibo stars over latitude and 
outflow velocity. There is a difference in the concentration toward 
the plane at low latitudes between stars with outflow velocities 
higher and lower than 10 kms -1 , respectively. Below 10 kms -1 
there is also a clear excess of sources (50%)in the Arecibo sam- 
ple with respect to the other samples in Fig. 4. These stars form 
possibly a thick disk. 



3.5 Radial and vertical variations 

Clearly, one double-exponential cannot describe the distri- 
bution of OH/IR stars. The older disk stars (few Gyr) seem 
to trace an exponential disk with rather large scalelength 
plus a central component with similar vertical scale (Table 
1). The younger stars (< lGyr) follow a more capricious pat- 
tern and have varying, but small vertical scales. This is not 
unexpected since they have completed only a few rotations 
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Figure 7 As Fig.4, for the Bulge AOSP sample ( "B" , 250 stars 
in the plot), the Disk AOSP sample ("1",165 stars) and the stars 
of the Arecibo sample with \b\ < 3° ("2P", 160 stars). There is 
still an excess of sources below V oxp =10kms- 1 for "2P". Note 
that the distribution of higher outflow velocities of "2P" is much 
more like "B" than like "1". 

around the galactic Centre since their formation and there- 
fore phase mixing, let alone radial mixing, has not yet been 
effective. There seem to be no sudden transitions in the dis- 
tribution; except that from the flat distribution (/ir — > cm) 
to the outer disk (/ir =1.5 kpc) for the high-outflow sources 
that is rather abrupt. As we will see later (Sect. 5), this may 
be around the radius of corotation with the pattern speed of 
the central Bar. The radial scale derived (via the same pro- 
cedure described earlier) from the Arecibo sample (using 
only |6| > 8° because of it's incompleteness in the plane) 
and from a sample of carbon-rich AGB stars (\£\ ~ 60°, 
\b < 15°, Loup et al. 1993), is ~ 10 kpc; the same as from 
the low-outflow OH/IR stars. So, except for the massive 
oxygen-rich objects, the scalelength is found to be similar 
for all AGB stars. The very different radial distributions of 
C-rich and O-rich AGB stars, respectively, is governed by 
the metallicity gradient only. 
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Outside the AOSP-survey window, we also derived 
scaleheights from the Arecibo subsample and the Loup sam- 
ple, and found h z ~ 500 pc for both. The full sample of 
outer-Galaxy AGB (not only OH/IR) stars of Blommaert 
et al. (1993) has h z Ss400 pc. The vertical scale appears to in- 
crease with latitude from ~100 pc to ~500 pc, irrespective 
of the longitude range. The increase in scaleheight within 
the Disk AOSP sample agrees with the diffusion models by 
Wielen (1977) for age increasing from ^1 Gyr to ^2 Gyr. 
The scaleheight of 500 pc for the oldest AGB stars is the 
same as that for white dwarfs as given in Mihalas & Binney 
(1981). 

3.6 In short 

In summary, the AGB stars are distributed in the thin (old) 
disk with a scaleheight of 100 pc for the youngest AGB stars 
(Ssl Gyr). The scaleheight increases continuously to 500 pc 
for AGB of ^5 Gyr. The scaleheights for the Disk and the 
Bulge are the same. The carbon-rich and oxygen-rich AGB 
stars form one population; the differences between the dis- 
tributions of the two groups are governed mainly by the 
metallicity gradient. There seems to be a small population 
that has a distinctly large scaleheight. 

All trends in the AOSP sample can be seen in Fig.l in a 
pictorial fashion. The exact numbers for the scales in Table 
1 should be used with care, since our deprojection method is 
indirect. Nevertheless, from modelling distribution functions 
for various of the samples used here the same trends numbers 
emerge (Sevenster 1997). Also, the scales in the Bulge region 
agree very well with those recently found from DENIS data 
(Gilmore priv.comm.). 

4 THE BAR 

4.1 Surface density 

In this section, we compare parametrized bar models to the 
surface density of the AOSP sample. We do not optimize fits 
quantitatively, since the choice of a parametrized bar model 
is already an arbitrary one. Rather, we take bar models that 
are known to give good approximations to other observations 
and see how well they agree with ours. 

The viewing angle <f> is the angle between the Bar's ma- 
jor axis and the line of sight toward the galactic Centre 
(Fig. 8(a)). In Fig.8(b,c) we show the surface density distri- 
bution as a function of longitude, N(£), for two fiat (two- 
dimensional) elliptical bars with gaussian density distribu- 
tion (cf. G2 model in Dwek et al. 1995). The viewing an- 
gles are taken to be 20°, as suggested by some observations 
(Dwek et al. 1995). The form of N(£) depends upon the dis- 
tance out to which we integrate or, in observational terms, 
the distance d out to which the data sample the Galaxy, so 
N(l) = N(£;d). For values of d ~R , N(l;d) essentially 
looks like the distribution arising from an rn = 1 distor- 
tion (lopsided; see Fig.8(a)), with its maximum toward pos- 
itive longitudes. With d = oo, however, the distribution is 
skewed toward negative longitudes. This is the result of the 
line of sight through the m = 2 distortion being longer on 
the far side than on the near side, for small values of abso- 
lute longitudes (Fig. 8(a)). This effect was first predicted by 
Blitz & Spergel (1991) and is also seen in the micro-lensing 



maps calculated by Evans (1994). The strength of this ef- 
fect obviously depends upon the parameters of the density 
distribution of the bar. 

We will try to find this effect in our Bulge AOSP 
sample. The double-peaked OH/IR stars are divided into 
two equally-sized samples by outflow velocity, at V" e x P =14 
kms -1 , with average V cxp of 11.3 (sample I) respectively 
18.3 (sample II) kms -1 (cf. Fig.l(e,d)). This gives a factor 
of 1.7 difference in stellar luminosity L», even if we assume 
fii = 2/j.u (equation (1)). Blommaert et al. (1997) find a 
range in /i of ~ 2 in the GC with IR observations. The range 
of masses in the Bulge (Sect. 2) makes a factor of 1.7 differ- 
ence in L* (~1.5 in Mi) very well possible. Since the OH 
masers are saturated, the OH luminosity Loh increases, on 
average, linearly with L*. The lower limit to the flux density 
5oh is the same for both samples, so the average limiting 
distance of sample II is a factor ~ 1.3 larger than of sample 
I. We can thus use the two samples to mimic the different 
integration limits in Fig. 8 (see also Sect. 2). 

The effect of skewed distributions should be clearest in 
the inner regions of the Galaxy (see Fig.8(b,c)). The ratios 
of the number of stars with 0° < I < 4° to the number of 
stars with 0° > £ > -4° are 39/35 (sample I) and 22/35 
(sample II). These ratios are in accordance with the theo- 
retical results shown in Fig. 8. To define these trends in a 
more sophisticated way, we sorted both samples on absolute 
longitude and calculated the cumulative sums of the sign of 
the longitude X/WI^!)' we ac ^ or subtract 1 for each star 
(Fig. 9). An axisymmetric distribution gives a line that hov- 
ers around zero. If negative (positive) longitudes are 'over- 
populated' the sum will steadily decay (rise). This sum is 
shown in Fig. 9 for the two data sets and for the bar model 
shown in Fig. 8(c), as well as for an N-body model (Fux 
1997) found to represent the AOSP sample well (Sevenster 
et al. 1999). The dotted curves give the 95% confidence limits 
for deviation from axisymmetry. Sample I (solid curve) never 
deviates significantly from axisymmetry, although there are 
local trends similar to those of the models with the inter- 
mediate cut-off. Sample II (dashed curve) , however, lies at 
or outside the 5% confidence level and coincides remarkably 
well with the models without distance cut-off. The set of 
stars with low V eX p would have an average distance cut-off 
of 9-9.4 kpc according to these models, the set with high 
Kxp of around 12 kpc, beyond which there is no significant 
contribution from the bar to the integrated density anymore. 
This agrees very well with the difference of a factor 1.3 in 
average distance derived from the relation between Kx P and 
stellar luminosity. 

We estimate the disk contamination to be ^20% within 
5° of longitude from the GC (Sect. 2; Sevenster 1997). If an 
axisymmetric component contributes significantly in the in- 
ner degrees, the evidence for the existence of the Bar would 
only become stronger. The XX^/KI) distribution would not 
change with the subtraction of a projected axisymmetric 
distribution of any relative density, but N(£ < £u m ) would 
become lower. This would make the probability of the devi- 
ation coming from a binomial distribution even smaller. In 
other words, the dotted curves in Fig. 9 would shift horizon- 
tally to the right, but the data-lines would remain in place 
on average. 

The evidence presented here is the first large-scale mor- 
phological evidence for a galactic Bar that cannot also be 
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Figure 8 a (left) Schematic view of the Bar from the galactic North pole. The viewing angle <j> is indicated; the near end of the Bar 
is most certainly at positive longitudes. The lines of sight (si and s2) at equal but opposite longitudes differ in length. The dotted 
arcs show schematically possible observation limits. The shaded circle represents a lopsided distortion, that would project to a similar 
surface-density distribution as the Bar for an observation limit of Kq (dotted arc 1). In b and c we show the total number of stars 
(arbitrarily normalized) along the line of sight in a two-dimensional elliptical barmodcl with gaussian density distribution for three 
different integration limits (8 kpc (1), 9 kpc (2) and oo (3); Rq =8 kpc). The integration limits correspond to the observational limits 
shown in a. For reference a gaussian is shown (dashed). Note the offset of the maximum toward negative longitudes for the largest 
integration limit, b (middle) Axis ratio 5=0.6, semi-major axis a=3.5 kpc (=3<T ga us). c (right) Axis ratio q=0.4, semi-major axis a=2.5 
kpc (=3o- gaus ). The viewing angle <j> is 20° in both panels b and c. The model in c is based on the K-band G2 model for i? m ax = 2.4 
kpc of Dwek et al. (1995). 
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Figure 9 a (left) The cumulative sum y](l/|l|) versus \t\ after sorting on \t\ for the two samples (thick solid and dashed lines) and for 
the G2— bar model from Fig. 8(c) with different integration cut-offs (short dashed lines, see Fig. 8). The dotted curves indicate the values 
for which the probability of the sum arising in a binomial distribution being larger (or smaller for negative values) than that value is 
5% (ie. single— sided). For \(\ < 0.5° and \l\ > 5° contributions by additional features influence the distribution (see Fig.l; Sect. 5). b 
(right) Same as a, but for a bar model consisting of two exponential elliptical— bar profiles with major— axis scalelengths of 200 pc and 
1 kpc (central— density contrast 8.5), with viewing angle of 45° and axis ratio of 0.5. The profile is an approximate fit to the density in 
the plane of an N-body bar and the viewing angle of 45° optimizes the fit between the AOSP data and the N-body model (Sevenster 
et al. 1999; Fux 1997). 



explained by a physical lopsided density distribution (Sev- 
enster 1996). The effect is also clearly visible in Fig.l(b,c) 
(sample 11,1) in a qualitative way. The same trend in the lon- 
gitude distribution is reported by Unavane & Gilmore (1998) 
who use L-band data in the plane. However, the necessary 
extinction corrections are estimates and of similar order as 



the asymmetries one is after. Unavane & Gilmore (1998) 
find that the Dwek models (E2, G3) give too large a devia- 
tion from axisymmetry, whereas the G2 model gives a good 
representation of our data (Fig.9(a)). The results in Fig. 9 
demonstrate that only for a given density model the param- 
eters of a bar can be optimized (Zhao 1997), as both rather 
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different models give good representations of the data. 
4.2 The Bulge's kinematic type 

The anisotropy parameter (Binney 1976) is a measure of 
the rotational support of a system. It is the ratio of the 
maximum rotational velocity of a bulge (U m ) to its central 
velocity dispersion (cr ). In combination with the flatten- 
ing e, it provides a way to distinguish between rotationally- 
supported, dispersion-supported and streaming-dominated 
systems (Illingworth 1977). The maximum rotation velocity 
of the Bulge occurs at £ ~ —18°, where the mean-velocity 
curve levels out (Fig.10), so V m = 140±20kms~ 1 . For cro, 
we use 153 kms" 1 (Blum 1996), which is more reliable than 
our determination of this second-order moment (135±20 
kms' 1 , Fig.lO(b)). This gives a value of V m /a of 0.9±0.1 
. Together with an ellipticity of e ~ 0.4 (Table 1; Dwek 
et al. 1995, GO- model; Kent 1992), this locates the galac- 
tic Bar in the Vm/oo _ e diagram between the oblate and 
the SAB bulges (Kormendy 1993), governed by somewhat - 
more-than-rotational support. The Galaxy fits well in the 
relation as an average SAB galaxy. 

Earlier values for the Bulge's dispersion, from measure- 
ments toward Baade's window, were much lower (113^5 
kms -1 ; Sharpies, Walker & Cropper 1990). Adopting this 
as cro , the Bulge would be located in the region of extremely 
triaxial bulges in the Vm/cro - e diagram. Baade's window 
(b — —4°) is too far from the plane to sample the central 
Bar dispersion. We conclude that the Bar is so flat that it 
cannot be assumed to be the dominant contributor to the 
distribution along the line of sight toward Baade's window 
(Table 1; see also Sevenster et al. 1999). Paradoxically, us- 
ing the dispersion in Baade's window one would obtain an 
anisotropy parameter indicative of a strong bar. 

5 RESONANT STRUCTURES 

In this section, we will concentrate on local features in the 
distributions of various tracers. For the AOSP sample, these 
structures are seen in Fig.l, marked as R_i, R-2 and R-3. 
We complement our own evolved-stellar data set with 2.4- 
GHz-continuum observations (Duncan et al. 1995) and a 
sample of star-forming regions (Comeron & Torra 1996). 
The synchrotron radiation dominating the continuum at 2.4 
GHz, is found to be a good tracer of the current locations 
of density waves (eg. Tilanus & Allen 1989). 

Duncan et al. (1995) present their data separated in 
a small-scale- (sub-degree) and large-scale distribution 
(Fig. 11). The small-scale emission arises largely in super- 
novae; the large-scale emission traces the molecular gas and 
density waves. For the uninitiated, an introduction to non- 
axisymmetric galactic dynamics can be found in Binney & 
Tremaine (Ch.6), and an introduction to terminolgy and lit- 
erature in Sevenster (1997). 

5.1 The corotation region 

The only direct method to find corotation (CR), via the 
pattern speed, is in most cases inapplicable (see Tremaine 

6 Weinberg 1984; Merrifield & Kuijken 1995). However, 
a number of indirect indicators can be used. A signa- 
ture of CR is often observed in the form of dust lanes, 



spiral-arm bifurcations or local decreases in the density (eg. 
Elmegreen 1996). A flat part of the rotation curve may out- 
line roughly the region between the inner-ultra-harmonic 
resonance (IUHR) and CR (Wozniak & Pfenniger 1997). 
This is observed in the edge-on SO-galaxy NGC 4570 (van 
den Bosch & Emsellem 1998). Similarly, the orbital stochas- 
ticity in this region (Contopoulos & Grosb0l 1986) will cause 
efficient radial mixing, smearing out the density gradient 
(see also Fux 1997). An easily observed result of the dynam- 
ics in this region are "inner rings" that are prominent in the 
majority of barred galaxies (see Buta 1996). 

All these features are in fact present in the Galaxy. The 
narrow gap between the so-called 3-kpc and Norma arms 
(Fig. 11), also seen in other tracers (eg. CO, Bronfman 1992), 
could well indicate CR. In good agreement with this is that 
the flat part of the rotation curve found for the AOSP sam- 
ple (Fig.10) stretches from I ~ -18° to £ ~ -25°. Those 
longitudes would, according to Wozniak & Pfenniger (1997), 
indicate IUHR and CR, respectively. The very flat density 
found for this region in Sect. 3 further reinforces this picture 
(Contopoulos & Grosb0l 1986). It also rules out the pos- 
sibility that the rotation curve is flat due to a logarithmic 
potential (density oc RT 2 ). 

Most importantly, we argue that the 3-kpc arm is the 
projection not of a spiral arm but of an inner ring, such 
as mentioned above, for the following reason. A subsample 
of the AOSP sample follows exactly the longitude-velocity 
structure of the 3-kpc filament (Fig. 12). A group of nine 
young (high-outflow) stars, that trace the kinematic struc- 
ture between 0° > I > —10°, stands out in the left panel of 
Fig. 12. These nine stars are at very low latitudes and have 
a very high median outflow velocity of 17.5 kms^ 1 . If the 
metallicity is the same as in the rest of the disk, the initial 
masses of this 3-kpc sample would be ~ 6 Mq and the ages 
~ 600 Myr (Appendix A). As those stars remain close to the 
gas for several galactic years, their trajectories must follow 
closely that of the gas, In other words, the gas filament must 
outline closed orbits instead of a (temporary) spiral density- 
wave maximum. With proper motions for these stars one 
could constrain the motion along the 3-kpc arm completely, 
something that is impossible to achieve with gas. 

In Fig. 13(a), we show longitude-velocity trajectories 
constructed by Mulder & Liem (1986). They solved gas- 
dynamical equations in a weakly-barred potential, using 
a finite-difference, hydrodynamic grid code. The gas is 
described by the inviscid Euler equations. The curves in 
Fig. 13(a) represent their preferred model for the 3-kpc arm, 
the main feature used to constrain the scales of their mod- 
els. The viewing angle is <f> =20° and Rcr = 0.98 Rq. In 
Fig. 13(b), the model spiral is seen face-on. Although this 
spiral reproduces the 3-kpc filament around £=0° very well, 
the tangent point to the 3-kpc arm is at much too high longi- 
tude, £ = -27° (Fig.l3(a)). Mulder & Liem (1986) show the 
longitude-velocity diagram for another model, for <fi =40° 
and -Rcr = 0.56 Rq. In this diagram the tangent point is 
much closer to the observed point (see Fig. 12) and also the 
central CO-disk's velocity signature is matched much bet- 
ter. The trajectories of both ILR arms remain similar in the 
inner regions. Unfortunately they do not give trajectories for 
this model but they note that the arms are "more concen- 
trated around CR". This would be in agreement with our 
inner ring, discussed above. The Mulder-Liem models are 
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Figure 10 The mean line— of— sight velocity (dash) in the inertial frame (for Vlsr =200 kms -1 ) smoothed data (double— peaked stars 
only) at b = 0° (a) and b = 2° (b). The gridsize is 1° by 6 kms -1 , the initial kernel sizes are 1° and 30kms _1 , respectively. We show 
errors derived via bootstrapping (Press et al. 1992, sec Sevenster 1997) for the features that we discuss in the text. (The mean velocity 
is the 50% value and dispersion is half the difference of the 83% and the 17% values of the smoothed velocity profiles.) 
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Figure 11 The observed surface-density distribution at 2.4 GHz due to Duncan et al. (1995). The features R_i, R_2 and R_3 are 
indicated (see Fig.l) as well as the "Norma arm" (tangent point —32°) and the "3-kpc arm" (tangent point —22°). For the large-scale 
(> 1°) distribution, the ten dashed contours are spaced linearly between 10% and 100% of the maximum (large-scale) density. For the 
small-scale (< 1°) distribution, the ten solid contours are spaced logarithmically between 3% and 95% of the maximum (small-scale) 
density. The large-scale emission is a good tracer of density waves. 



forced to be quasi steady. Therefore, gas flows around CR 
would not be represented correctly and it could well be that 
the two arms in Fig. 13(b) connect in reality to form a ring. 

Binney et al. (1997) find, deprojecting the COBE map 
with imposed eight-fold symmetry, two density enhance- 
ments on the minor axis at 3 kpc from the Centre (£ — 
+17°, —22°). They suggest these are the L^s points. With 
four-fold symmetry, however, they find a, leading, spiral. 
The longitude £ = —22° of their feature coincides with the 
maxima in the OH/1R stars and the 2.4 GHz emission used 
in this section. As Binney et al. (1997) note, there may be a 
significant contribution of young stars to the K-band surface 
density (Rhoads 1996). We argue that their density features 
are incorrect deprojections of the inner ring. In any case, the 
derived loci for CR are very similar. 

The maxima R_2 and R_3 interestingly border those 
of the 2.4 GHz emission (Fig.l). Also, the locations of the 
maxima in the older and younger OH/IR stars, respectively, 
are slightly different. Such displacements may be caused by 



the intricate streaming and diffusion processes around CR 
(Roberts, Huntley & van Albada 1979; Kenney & Lord 1991; 
Vogel et al. 1993; Tilanus & Allen 1989), but it is virtually 
impossible to draw conclusions from them. 

5.2 The inner— Lindblad resonance 

Other than a corotation resonance and outer-Lindblad reso- 
nances, the existence of an inner-Lindblad resonance (ILR) 
in a galaxy is dependent on its exact potential. If it exists, 
is may show via a a variety of indicators. The ILR is usually 
outlined by a "nuclear ring", often accompanied by massive 
star formation (see eg. Buta 1996; Phillips 1996). A so-called 
double-wave feature in the stellar rotation curve (Bettoni 
1989) may exist as a result of orbits trapped around the 
the retrograde x± orbit family inside ILR (Pfenniger 1984; 
Wozniak & Pfenniger 1997; Contopoulos & Papayannopou- 
los 1980) . Finally, if an ILR exists, the gas flowing toward the 
centre, following bar-induced instabilities, will follow offset 
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Figure 12 The longitude-velocity diagram for HI (grey scale) in the plane and AOSP stars with \b\ < 1.5°. Ovcrplotted in the left 
panel are stars with V e xp> 14 kms -1 (the young sample) and in the right panel stars with V CX p< 14 kms" 1 (the old sample plus the 
single— peaked stars). The 3-kpc arm is clearly seen in the HI as the dark linear feature extending from (+5°, —20 kms - 1 ) to (—20°, 
-130 kms - ). The young stars have a subgroup of nine massive stars that follows exactly the longitude— velocity structure of this arm. 
The old sample shows no connection at all, and deviates in general much more from the HI longitude-velocity diagram. (Courtesy A. 
Kalnajs) 



spiral arms (Athanassoula 1992). In absence of an ILR, the 
gas flows to the centre directly along the major axis of the 
bar. 

Comeron & Torra (1996) report an elongated ring- 
shaped maximum in the deprojected density distribution of 
ultra-compact HII regions (UCHII). They tentatively iden- 
tify this with an ILR- or nuclear ring and give radii of 1.3 
kpc and 1.9 kpc, for the southern- and northern hemisphere, 
respectively (corrected for R =8 kpc). In general, these nu- 
clear rings are found to have radii around 1—1.5 kpc (Buta 
1996; Freeman 1996). 

The radio-continuum distribution (Fig. 11) does not 
show an increase in the number of small-scale sources at 
I ~ — 9? 5 (1.3 kpc). In that case, star formation must have 
started too recently for massive stars to become supernovae 
(~ 10 7 yr) and star formation must have ignited between 
10 7 yr and 5xl0 4 yr ago (the ages of the UCHIIs). 

A double-wave feature in the OH/IR-star rotation 
curve can indeed tentatively be identified in Fig. 10(a) (I ~ 
+3°,— 4°). As the AOSP sample is the only large-scale 
stellar-kinematical sample at low latitudes, it is impossible 
to verify this feature with other data. According to Woz- 



niak & Pfenniger (1997), the deepest minimum in the ro- 
tation curve arises around 0.25-Rcr, when the retrograde 
orbits remain concentrated around the X4, family. It is not 
straightforward to say where the maximum effect will take 
place in projection; a simple tangent-point assumption is 
probably not valid. The Max(i5)-model of Wozniak & Pfen- 
niger (1997, see their Figs. 4, 5) is a good example. In that 
model, the maximum effect, from the viewing point of the 
Sun, of this "counter-streaming" comes from the minor axis 
of the bar. A location on the minor axis at 0.25-Rcr. would 
indeed be seen roughly at the longitudes of the observed 
double-wave feature. 

In external barred galaxies, this feature is only seen 
at intermediate inclinations (Bettoni & Galletta 1997); the 
Galaxy might be the only edge-on system that allows detec- 
tion of the double wave, as the disk-foreground contamina- 
tion is relatively small with respect to other galaxies. Some 
early-type galaxies are known to harbour kinematically- 
distinct cores (see eg. Carollo et al. 1997) that can give a 
similar feature in the rotation curve. The cores contribute 
a much larger fraction to the total galaxy than do the ret- 
rograde orbits mentioned above and so the feature is not so 
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Figure 13 a The longitude-velocity diagram for AOSP stars with \b\ <1° (open circles). The filled circles are the 7 stars from R_2 with 
\b\ < 1° and the filled triangles the 9 stars connected to the 3-kpc arm. The solid curve is the "near ILR arm" from a model by Mulder 
& Liem (1986) with <j>=20° and i?cjj=0.98 Rq, the dashed curve the "far ILR arm" (see b). The velocities are in the inertial frame 
(changed from Mulder & Liem (1986) for Vl SR =200 kins" 1 instead of 250 kms -1 . b The ILR arms from Mulder & Liem viewed from 
the galactic North pole. The thickened parts indicate roughly the positions of the two groups of stars marked in a if they are connected 
to those spiral arms. Note that the tangent point to the 3-kpc arm is at too high longitude for the model arm. 



easily masked by the foreground disk. We do not think the 
feature observed in the Galaxy is caused by such a decoupled 
core because the wave's velocity maximum and minimum are 
small with respect to the rotation at large radii. We checked 
that the feature is not the result of smearing out the ef- 
fect of the rapidly-rotating galactic-centre OH/IR stars (see 
Sect. 3. 2; Lindqvist et al. 1992; Sjouwerman et al. 1998a). 

An offset gas flow around ILR would be virtually impos- 
sible to detect in the edge-on jungle of the galactic disk. In 
Fig. 13(a), next to the 3-kpc stars, another group of stars is 
highlighted. The 13 stars forming the feature R_i (Fig. 1(e)), 
seen in the low-outflow, older OH/IR stars only, follow a co- 
herent filament in the longitude-velocity diagram that co- 
incides with part of the far ILR arm (see Fig. 13(b)). All 
13 stars lie on this filament, although in Fig. 13 we have 
plotted only 7 with \b\ < 0.6°. (Note that these stars were 
selected before overplotting them on the ILR-arm model.) 
In Fig. 13(b) we indicate roughly where the two small groups 
of OH/IR stars would be located in the Mulder-Liem model 
(see Sect. 5.1 for caveats concerning this model). The com- 
pactness (see Fig. 1(e)) of this group of relatively old stars, 
can indeed be understood more easily if they are spread out 
along the line of sight as shown. However, the far-advanced 
age of this group of ILR-arm OH/IR stars is puzzling and 
prevents us from accepting their location in this picture too 
readily. 

Exemplary of the ILR region in our Galaxy may be 
the inner region of the galaxy IC4214 (Hubble type SABab, 
Saraiva 1996). 



6 THE PRESENT SHAPE OF THE INNER 
GALAXY 

The evidence presented in the previous sections may be com- 
bined to give the following description of the inner Galaxy 
(Fig. 14). There is a stellar Bar (Sect. 4.1) that is fairly flat 
and weak (Sect.3;4.2). It corotates at ~3.5 kpc and is sur- 
rounded by an inner ring, inside corotation, roughly between 
2.5 kpc and 3.5 kpc (Sect. 5.1). Possibly an inner-Lindblad 
resonance exists, with corresponding nuclear ring, between 1 
kpc and 2 kpc (Sect. 5. 2). We suggest that outside corotation 
a four-armed, trailing spiral exists, which is a bifurcation of 
a two-armed spiral mode. As such, the Galaxy would be sim- 
ilar in appearance to NGC 1433 (type (Ri)SB(rs)ab; Buta 
1995; for a picture see Buta 1996). 

If we assume the rotation curve of the axisymmet- 
ric part of the Galaxy is described by 184 km s" 1 kpc -1 
R 01 (Allen, Hyland & Jones 1983), a pattern speed of 60 
kms -1 kpc -1 is implied for 7?cr=3.5 kpc. The radii of ILR, 
IUHR and OLR in this case are 0.85 kpc, 2.2 kpc and 6.3 
kpc, respectively. The radius of ILR depends crucially on the 
adopted potential and larger radii are easily possible in more 
centrally concentrated potentials. The radius of IUHR, for 
60 kms -1 kpc -1 , agrees well with our suggestion that the in- 
ner ring, discussed in Sect. 5.1, lies between IUHR and CR. A 
sudden decrease in metallicity is observed around R — 6 kpc 
(Simpson et al. 1995). If this marks the end of the zone-of- 
infiuence of the Bar's mixing, as suggested by the authors, 
one expects this to be around OLR. This agrees with the 
radius of OLR for a pattern speed of 60 kms -1 kpc -1 . 

Comparison with published values for the corotation 
radius and pattern speed proves little conclusive. Values for 
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Figure 14 Top-view sketch of the Galaxy as described in Sect. 6. 
The larger ellipse outlines roughly the bar density; the smaller the 
ILR. The two open circles inside the ILR indicate the loci where 
the double-wave rotation feature may originate. The short curves 
emerging from ILR are part of the ILR arms (Mulder & Liem 
1986) connecting to IUHR— CR. The thin, solid circle around the 
bar is CR. The thick dashed ring indicates roughly the IUHR, 
the inner ring lies between the IUHR- and the CR circle. The 
spiral arms outside CR are inspired by the picture by Georgelin 
& Georgelin (1976), structured by the assumption that the four 
spiral arms are two pairs, bifurcating from one two-armed spiral 
inside CR. The black dot marks the position of the Sun; the figure 
is approximately to scale with the distance of the Sun to the 
galactic Centre (8 kpc). Compare to Fig. 6 in de Vaucouleurs & 
Pence (1978), who sketch a similar view of the Galaxy based on 
their derived Hubble type (SAB(rs)bc). 

corotation range from 2.4 kpc (Binney et al. 1991), 3.5 kpc 
(Weiner & Sellwood 1996; Englmaier & Gerhard 1999), 4.8 
kpc (Fux 1997) to 9 kpc (Amaral & Lepine 1997), and, sim- 
ilarly, for the pattern speed from 19kms _1 kpc -1 (Wada et 
al. 1996) to 118 kms^kpc" 1 (Yuan 1984). 

6.1 Elongation 

In Sect. 4. 2, we found that the kinematic type of the Bar in- 
dicates a rather weakly-barred distribution. This is in agree- 
ment with the large in-plane axis ratio of the density, found 
in several studies : q p = q$,n = 0.6 (Binney et al. 1997, stel- 
lar density) or q p = 0.7 (Wada et al. 1996, gas kinematics), 
although Nikolaev & Weinberg (1997) find q p — 0.4. From 
Fig. 9, q p ~ 0.5 is not excluded by our data. Due to mix- 
ing (Friedli, Benz & Kennicutt 1994), a correlation between 
the in-plane axis ratio of bars and the metallicity gradi- 
ent of their host galaxies exists. From an empirical relation 
between bar axis ratio and [O/H] gradient (Martin & Roy 
1994), one finds q p = 0.7 for a gradient of —0.07 dex kpc -1 
(Smartt & Rolleston 1997). 

Inner rings are usually elongated along the bar and fol- 
low the shape of the main closed orbits; their axis ratio is 
q^i (see Buta 1996) and thus - in general - larger than 
q p . The elongation can be constrained by the radial motion 



along the ring. At £ —0°, the line-of-sight velocity is —53 
km s _1 , or in other words a radial velocity Vr of +53 km s _1 
(Fig. 12; Oort 1977). If this is the line-of-sight component 
of a streaming along the ring, Vji, we can define an angle 
/3 = acos(Vk/Vj|). If we assume further that the ring is an 
ellipse, with a viewing angle (j> of the major axis (near end at 
positive longitudes), and that the streaming is along the el- 
lipse, we can derive the axis ratio q e via ql = tan(/3— 0) tan 4>, 
with < 4>< (3 < f . For Vj| = 200 kms" 1 , /3 = 75° and we 
find q c = 0.76 for <j> = 45° (q c = 0.75 for </> = 25°; q c = 0.61 
for 4> = 65°). For Vj| = 150 - 300 kms" 1 , (3 = 70°-80° and 
q c = 0.68 - 0.84. Hence, for 20° < </> < 55°, the observed Vr 
is consistent with q e = q&i > 0.6 for a wide range of Vy. The 
largest q& arises for <f> = 0.5 /3. 

In general q<s> > g< 6 / > q p , so the radial velocity of the 
3-kpc arm can be accommodated in an almost-round poten- 
tial (see also Sil'chenko et al. 1997), contrary to some claims 
(eg. Mulder & Liem 1986). 

7 EVOLUTION OF THE INNER GALAXY 

The morphological properties of bars are found to be re- 
lated to their formation mechanism (Noguchi 1996), to the 
rotation curve (Combes & Elmegreen 1993) and to the Hub- 
ble type of the host galaxy (Elmegreen & Elmegreen 1985). 
Early-type-galaxy bars (SsSb in this context) have flat den- 
sity profiles and an ILR, end close to their corotation radius 
and form via strong tidal interaction with other galaxies. 
Late types (<; Sc) have exponential density profiles, no ILR, 
end well inside (sometimes halfway) their corotation radius 
and form via disk instabilities. 

In this section, we will speculate on the formation and 
evolution of the galactic Bar, and on how it influenced the 
structure of the inner Galaxy. 

7.1 How did the Bar form ? 

The Galaxy is usually classified as an Sb or Sbc galaxy in 
the Hubble sequence (de Vaucouleurs & Pence 1978) and is 
as such just intermediate to the two classes described above. 
The Bar's steep density profile (Sect. 3, 4; Binney et al. 1997) 
and length (Ren/a > Rcr/Rivhr = 1.4) indicate a late 
type. Is it possible that the Bar in our Galaxy indeed formed 
via disk instability ? Studies of Toomre's instability param- 
eter Q (Toomre 1964) show that the Galaxy is (marginally) 
stable at all radii (Fuchs & von Linden 1997; Lewis & Free- 
man 1989). Earn & Lynden-Bell (1996) also find that, for 
the Bulge's present rotation curve V c oc R 01 (p oc R~ 18 , 
Allen et al. 1983) , m = 2 modes will not be sustained via 
orbit cooperation (alignment). 

However, in the first gigayears of the life of the Galaxy 
(3>5 Gyr ago), when the stellar dispersions had not yet in- 
creased to their maximum value (Freeman 1991), the disk 
was probably unstable (Fuchs & von Linden 1997). Espe- 
cially if there was a weak trigger (the Large Magellanic 
Cloud, Weinberg 1996) it is possible that the Bar formed via 
a disk instability, more than 5 Gyr ago. The strong correla- 
tion between kinematics and metallicity (eg. Zhao, Spergel & 
Rich 1994) favours secular bar formation and thus excludes 
strong interaction as a mechanism, pointing to a disk-bar 
connection rather than to a halo-bar connection. The similar 
scaleheights for the two components found in Sect. 3 (Table 
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1) support a disk-bar connection and the formation of the 
bar from the disk via a planar instability. 

7.2 When did the Bar form ? 

In general, it is suggested that the Bar is a young component 
with respect to the spheroidal components (Bulge, Halo). 
Several authors find ages for the Bar of 6 Gyr to 9 Gyr 
(Ng et al. 1996; Gerhard & Binney 1993; Wyse, Gilmore & 
Franx 1997). The OH/IR stars in the direction of the Bar 
have an upper age limit of 7.5 Gyr (Sect. 2). If these OH/IR 
stars are the result of renewed star formation in response to 
the formation of the Bar (see discussion Sect. 7. 3), then the 
Bar must be at least 8 Gyr old. Such an age might allow 
for the disk-instability formation mechanism (see above). 
The scenario of Ng et al. (1996; their Fig. 13) explains the 
absence of very old Bulge OH/IR stars and the similarity of 
the metallicity of OH/IR stars in the Bar and in the Disk 
(Sect. 2). The absence of Bulge OH/IR stars older than 7.5 
Gyr indicates that the already existing old Bulge (>13 Gyr; 
Ng et al. 1996; Wyse et al. 1997) was not forming significant 
numbers of 1-6 M stars in the plane in this era. 

7.3 What happened next ? 

Of special interest is the bar-induced mass flow to the cen- 
tral regions of galaxies (Hasan, Pfenniger & Norman 1993; 
Wada & Habe 1992,1995) that has strong influence on the 
evolution of the bar itself, eventually resulting in its destruc- 
tion. Wada & Habe (1992, 1995) find that the presence of 
an ILR is essential, but however, Athanassoula (1992) and 
Piner, Stone & Teuben (1995) find that more generic mass 
transport always occurs, even preferentially in systems with- 
out ILR. 

The gas, flowing to the centre, causes instabilities and 
star formation. In order to induce Jeans' instabilities at the 
scales of giant molecular clouds (Mj < Mgmc, eg. Gerrit- 
sen & Icke 1997) 10 8 M of gas has to be transported to 
within 100 pc (for Mgmc = 1 x 10 6 M ; the masses of the 
SgrA,B GMCs are 2, 6 x 10 6 M , Stark et al. 1991). The ob- 
served gas-inflow rate into the central 100 pc is 0.01 M yr _1 
(von Linden, Duschl & Biermann 1993) and presumably was 
higher just after formation of the Bar (eg. Piner et al. 1995). 
This means it would take less than 10 Gyr, or more or less 
the estimated age of the Bar, for the central density to build 
up enough for the onset of central star formation. 

Evidence for massive, on-going star formation in the 
galactic Centre from eg. the presence of H2O masers has 
been negated (Sjouwerman & van Langevelde 1996). Never- 
theless, the rate inside 500 pc is higher than elsewhere, ~0.5 
Mq yr~ (Gusten 1989) and the OH/IR stars in the galactic 
Centre (inside 100 pc) appear to be younger on average than 
the global population (Sjouwerman et al. 1998b; Blommaert 
et al. 1997; Wood, Habing & McGregor 1998). 

The galactic-centre OH/IR stars have a bimodal distri- 
bution of outflows (Fig.4) that differs from that of the Bulge 
sample with more than 99.9% significance (Kolmogorov- 
Smirnov). The lower outflow velocities are distributed in 
a rather wide range (26% fractional scatter) around 11.4 
kms -1 , the higher in a narrow range (10%) around 19.4 
kms -1 . Especially when allowing for a scatter in the metal- 
licity, the masses of the high-outflow galactic-centre stars 



must span a relatively narrow range as well (Appendix A). 
This suggests a short episode of star formation, which ac- 
cording to Sjouwerman et al. (1998b) took place more than 
1 Gyr ago. The two groups of galactic-centre OH/IR stars 
are separated roughly by \4 xp = 15 kms -1 , corresponding 
to an age of approximately 1.5 Gyr (Sect. 2, Appendix A) for 
solar metallicity. The sudden horizontal step in the cumula- 
tive distribution around 15 kms" 1 suggests that around 1.5 
Gyr ago there was either a sudden increase of the ambient 
metallicity or of the average mass of the formed stars. Obser- 
vations show that the high-outflow galactic-centre OH/IR 
stars have metallicities of a factor of two higher than the 
low-outflow- and the Bulge objects (Blommaert et al. 1997; 
Wood et al. 1998). Hence, in the case of a burst starting 
1.5 Gyr before present, the progenitor masses of the present 
high-outflow population must have been well over 3 M 
(Bertelli et al. 1994) , in good agreement with Wood et 
al. (1998). For twice-solar metallicity, the luminosity of these 
objects would be ~90% of the average luminosity of the 
Disk AOSP stars (Bertelli et al. 1994) which, with the high 
outflow velocities, agrees well with equation (1). Sjouwer- 
man et al. (1998b) find that in total more than 10 7 M of 
stars would have formed in the starburst. However, they use 
half-solar metallicities in which case the luminosity would 
be ~110% of the Disk AOSP stars (still for ages of 1.5 Gyr). 
Such a model does not satisfy equation (1) at all. As this 
model uses lower progenitor masses, we conclude that in 
our high-metallicity scenario the burst would be a factor of 
five times more massive than estimated by Sjouwerman et 
al. (1998b). 

More recent (~ 7 Myr and ~100 Myr ago) and much 
less massive starbursts are reported by Krabbe et al. (1995). 
These are too recent to show up in the outflow-velocity dis- 
tribution of the OH/IR stars. 

7.4 The complete scenario 

The following speculative scenario emerges. The Bar formed 
~ 8 Gyr ago via a disk instability. Subsequent gas flows re- 
vived star formation in the Bulge region and the progenitors 
of the present-day Bulge OH/IR stars started to form. The 
metallicity in the inner plane became similar for all radii. 
Enriched material kept flowing into the central 100 pc, from 
the region inside corotation, and around 1.5 Gyr ago, a star 
burst ignited, leaving a significantly more metal-rich popu- 
lation in the galactic Centre. 

Another possible consequence of the build-up of mass 
in the centre is the formation of an ILR, typically when the 
mass inside 200 pc is 1-2% of the mass of the galaxy (Friedli 
& Benz 1993). The ring of UCHIIs (Sect. 5. 2) could indicate 
that an ILR exists in the Galaxy, and that has formed re- 
cently because it is not seen in any older (>10 Myr) popu- 
lation. The critical density for a nuclear-ring star burst is 
0.6k 2 /G (Elmegreen 1994a, b), yielding a required ring mass 
of SslxlO 8 M which could have accumulated from the in- 
flowing gas in ~1 Gyr, once the ILR formed (O.lMoyr -1 , 
Friedli & Benz 1993, for q p = 0.7). The ILR would hence 
have formed coevally with the central star burst derived from 
the OH/IR star population. 

The total mass within 100 pc is > 5 x 10 s M (Lindqvist 
et al. 1992; Kent 1992). This is ~3% of the total mass of 
the Bulge (~2x 10 10 M ; Zhao 1996; Blum 1995; Kent 
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1992). Such a concentration is sufficient to destroy the Bar 
completely (Hasan et al. 1993). 

The ILR may slow down the self-destruction of the Bar, 
by decreasing the rate of gas flow to the centre, but does not 
stop it (Hasan et al. 1993; Piner et al. 1995). In addition, in- 
side ILR no bar-supporting orbits would exist (Contopoulos 
& Papayannopoulos 1980). With already an almost-critical 
amount of mass in the central regions, the Bar may be com- 
pletely dissolved within a few Gyr. 

A view of the possible future of the Galaxy can be ob- 
tained by looking at NGC 7217 (Athanassoula 1996). This 
galaxy has an inner, an outer and a nuclear ring, but no clear 
bar. It may be a remnant of a once-barred galaxy that dis- 
solved, which could also account for the high fraction (25%) 
of counter-rotating orbits in this galaxy (Merrifield & Kui- 
jken, 1994). 

The observed CO parallelogram (Bally et al. 1988) was 
elegantly explained as the inner cusped x\ orbit by Binney et 
al. (1991). However, Liszt & Burton (1978) already showed 
that the parallelogram is most likely the result of a tilt in the 
inner CO distribution. In this case, the constraint, placed 
on the Bar's parameters by identifying this inner cusped 
xi orbit (yielding a very low viewing angle and corotation 
radius), can be lifted. In our scenario, the CO would be on 
X2 orbits inside ILR. 



8 CONCLUSIONS 

We have given several new arguments in favour of the tri- 
axiality of the central Galaxy. In the inner 10°, perspective 
effects provide evidence for the existence of a Bar. More indi- 
rect evidence comes from the presence of resonant features. 
We argue that an inner ring, with a flat density distribution 
and rotation curve, and a gap in the Norma-arm region in- 
dicate the region inside corotation. A possible nuclear ring 
and a double-wave feature in the stellar rotation curve are 
indicative of the existence of an inner-Lindblad resonance. 
We find that corotation is at 3.5 kpc, the bar ends within 
2.5 kpc (IUHR) and the pattern speed is ~ 60kms _1 kpc -1 . 
The radius of the ILR is most likely larger than what would 
be derived from a R 01 rotation curve (Allen et al. 1983). 

The much-studied 3-kpc arm is probably an inner ring, 
connected to the IUHR-CR region rather than OLR (Yuan 
1984; Binney et al. 1991) or ILR (Amaral & Lepine 1997). Its 
radial motion at £ =0° can be explained well by streaming 
along mildly elongated closed orbits (q > 0.7) for a wide 
range of viewing angles. 

The galactic Bar is no extraordinary specimen. The 
value of the anisotropy parameter Vm/oo and flattening e 
define it as an SAB-bulge. The Hubble type of the Galaxy 
is SAB(rs)b (and possibly an additional (nr) for the nuclear 
ring), based on its axis ratio, its spiral and ring structure 
and its kinematic type (cf. SAB(rs)bc, de Vaucouleurs & 
Pence 1978). 

Based on the maximum ages of the OH/IR stars in the 
direction of the Bulge, we speculate that the Bar's forma- 
tion took place ~ 8 Gyr ago. Induced mass flow to the cen- 
tre changed the mass distribution significantly. As a conse- 
quence, at ~ 1.5 Gyr ago a major star burst ignited (within 
100 pc) and an ILR formed. The Bar could be in the final 
stage of its existence. 



The OH/IR stars in our sample at low latitudes are all 
part of the intermediate-age Bar and do not trace the old, 
axisymmetric Bulge mentioned by eg. Ng et al. (1996) and 
Wyse et al. (1997). The Bar determines, via its influence 
on the metallicity gradient, the radial distribution of AGB 
stars in the disk (~0.5 to 1.5 R©), as the ratio of oxygen- 
rich to carbon-rich AGB stars is governed by the metallicity 
gradient. These two groups together form one population 
in terms of their galactic distribution. The AGB stars are 
distributed in the thin (old) disk with a scaleheight of 100 
pc for the youngest AGB stars ( ^1 Gyr) and 500 pc for 
AGB stars older than ^5 Gyr. There is possibly a disjunct 
1-kpc-scaleheight population of OH/IR stars, that is seen 
only very locally (I ~50°) and is between 1 Gyr and 10 Gyr 
old (depending on metallicity). 
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Appendix A The luminosity— metallicity— outflow- 
velocity relation 



In this appendix, we outline the derivation of the relation 
between (unknown) stellar luminosity, outflow velocity and 
envelope metallicity for OH/IR stars. 

We assume a star of luminosity L» and an isotropic and 
constant mass outflow M. At a distance Ro the gas tempera- 
ture has decreased to a value of T c , equal to the condensation 
temperature of some solid species. The particles that form 
are immediately accelerated to their final outflow velocity 
V a xp. Because of friction they take the gas along. We have 
now the following equation : 



MK, 



tL„/c 



(Al) 



as the momentum change of the outflowing mass equals the 
radiation pressure on the dust, r is the dust opacity of the 
circum-stellar envelope and c the speed of light. Also : 



M = AnRoPoVe, 



(A2) 



with po the density at radius Ro. From equation (A2) and 
the fact that Kx P is constant it follows that p — po(Ro / R) 2 ■ 
We then obtain : 



T = — ko 
P 



p(R)dR = -KopoRo 



Bo 



/' 



(A3) 



with ko the mean absorption coefficient and p the gas-to- 
dust ratio of the circum-stellar envelope. The dust forms 
when the radiation temperature has fallen to T c : 



aT* = L„/(4ttJ^) 



(A4) 



with a the Stefan-Boltzmann constant. We have four equa- 
tions and seven variables (M, L*,p, r, Ro, po, Kx P ) and thus 
three parameters can be chosen freely. We will take M, L, 
and p. It is then easy to show that : 



AL,p 



where the constant A is given by 



A = 



(4tt) 3 T c 4 c 2 



(A5) 



(A6) 



Thus V exp is independent of M. 

In a more elaborate model one may drop the assumption 
that the dust acquires its final velocity immediately. The gas 
is being dragged along by the dust; there will be a drift be- 
tween those components. Analytic solutions no longer exist 



and equation (A5) is slightly different (Habing et al. 1994) 
= AL*p~ 17 (A7 



(A7) 



By dividing a sample of OH/IR stars into two according 
to outflow velocity one thus effectively divides the sample 
into more and less luminous stars. All existing models of 
AGB-star evolution (eg. Vassiliades & Wood) predict that 
the more massive a main-sequence star is, the higher its 
AGB luminosity. Thus stars of high V cxp are younger that 
those of low V cxp . This is confirmed by the fact that stars 
with low 14x P have larger scaleheight than those with high 
K. xp (first discussed by Baud et al. 1981). 

In principle there could be a conspiracy. The stars with 
higher V cxp could have higher metallicities and lower val- 
ues of p. A variation in p could thus cancel the effect of a 
variation in V eX p and stars with higher V cxp would be older. 
However, for OH/IR stars, p tx Z" 1 (Habing et al. 1994). 
For the conspiracy to work, Z tx V e \p and Z must differ by 
at least a factor of four between stars to explain the observed 
outflow velocities (<10 to >20 kms -1 ). It is therefore un- 
likely that the conspiracy should influence our conclusion, 
because the OH/IR stars are younger than ~ 7 Gyr and 
oxygen-rich and thus have all formed in similar metallic- 
ity regimes (different by less than a factor of three over 7 
Gyr). A factor of three difference in Z may exist between the 
extremes of the outer-Galaxy and the metal-rich galactic- 
centre populations. 



Appendix B Turn-off from the AGB— evolutionary 
track 



In this appendix, we discuss a method to determine upper 
age limits for an ensemble of OH/IR stars from their "turn- 
over" [25]- [60] colours in the IRAS two-colour diagram. 

According to Garcia Lario (1991), OH/IR stars evolve 
along the evolutionary track (van der Veen & Habing 1988) 
in the IRAS two-colour diagram (Fig.Bl) until they reach 
maximum colours (R^jR^) dependent upon their main- 
sequence mass. At that point, they leave the evolutionary 
track to evolve with approximately constant R32 = R32 to- 
ward higher R21. The higher -R32, the higher the initial mass 
of a star : 



log Mi = 



R> 



2.42 



0.72 



- 2.45 /3.2 



(Bl) 



From Fig.Bl it is seen how one can thus determine the upper 
and lower limit to the initial masses and ages in a sample. 
The largest turn-over colour R^ 2 is indicated by B and the 
smallest by A. With equation (Bl) we find the highest and 
lowest masses present in the sample and with the isochrones 
by Bertelli et al. (1994) we find the corresponding lowest and 
highest ages. The lowest age thus found is an upper limit to 
the true lowest age as there may be younger stars that have 
not yet turned away from the evolutionary track (Fig.Bl). 



This paper has been produced using the Royal Astronomical 
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Figure Bl The IRAS two-colour diagram (van der Veen & 
Habing 1988). The colour R21 equals [12]-[25] which is defined as 
2.5 log S , 25 /J m/<S'i2 / jm and accordingly R32 = [25]-[60]. The solid 
curve indicates the OH/IR-star evolutionary track; the dashed 
lines outline the various regions as they were defined by van der 
Veen & Habing (1988). The two dotted lines give the largest and 
smallest R^ 2 (see text) for this fake sample. 



